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We have performed a search for the rib(lS) meson in the radiative decay of the T(2S) resonance 
using a sample of 91.6 x 10 6 T(2S) events recorded with the BABAR detector at the PEP-II B factory 
at the SLAC National Accelerator Laboratory. We observe a peak in the photon energy spectrum 
at E 7 = 609.31^ (stat) ± 1.9(syst) MeV, corresponding to an rj b (lS) mass of 9394.2^;® (stat) ± 
2.0(syst) MeV/c 2 . The branching fraction for the decay T(2S) — > jr/ b (lS) is determined to be [3.9± 



l.l(stat)l£;g(syst)] x 10 . We find the ratio of branching fractions B[Y(2S) jr) b {lS)]/B[Y(3S) 
7776(15)] = 0.82 ± 0.24(stat)t£ig(syst). 

PACS numbers: 13.20.Gd, 14.40.Gx, 14.65. Fy 



A candidate for the 775(15) meson, the ground state of the bottomonium system, was recently observed in the 
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radiative decays of the T(3S) The BaBAR experi- 
ment has accumulated a large sample of data at the peak 
of the T(2S) resonance, where radiative T(2S) decays 
are also expected to produce the %(15) meson. Theo- 
retical predictions for B[T(2S) — ► jrn,(lS)] range from 
(1 - 15) x 10" 4 @. A 90% confidence level upper limit 
of B[T(2S) 7??h(15)l< 5.1 x 10~ 4 is provided by the 
CLEO III experiment [H]. 

The ratio of branching fractions for the transitions 
Y(2S) -> 7776(15) and Y(3S) -> 77? 6 (l£) is dependent 
upon the overlap integrals of the relevant bottomonium 
wave functions 0, enabling a test that the observed 
state is the r]b(lS) meson. More generally, the mea- 
sured hypcrfine mass splitting between the triplet and 
singlet states in the bottomonium system provides a bet- 
ter understanding of nonrelativistic bound states in QCD 
and the role of spin-spin interactions in quarkonium mod- 
els di. 

In this Letter, we report evidence for the radiative 
transition T(2S) — > jj]b(lS). Hereafter r]b(lS) will be 
abbreviated as 775. 

The data used in this analysis were collected with 
the BaBAR detector Q at the PEP-II asymmetric-energy 
e + e~ storage rings. The primary data sample consists of 
14 fb _1 of integrated luminosity collected at the peak of 
the Y(2S) resonance. An additional sample of 44 fb _1 
collected 40 MeV below the T(iS) resonance is used for 
background and efficiency studies. The trajectories of 
charged particles are reconstructed with a combination of 
five layers of double-sided silicon strip detectors and a 40- 
layer drift chamber, both operated in the 1.5-T magnetic 
field of a superconducting solenoid. Photons are detected 
with a CsI(Tl) electromagnetic calorimeter (EMC). The 
photon energy resolution varies from 3.4% (at 300 MeV) 
to 2.8% (at 800 MeV). Hereafter we quote values of E 7 
measured in the center-of-mass (cm.) frame. 

The signal for T(2S) — > 7 77k is extracted from a fit to 
the inclusive photon energy spectrum. The monochro- 
matic photon from this decay should appear as a peak 
in the photon energy spectrum near 615 MeV on top 
of a smooth nonpeaking background from continuum 
(e + e~ — * qq with q — u,d, s, c) events and bottomonium 
decays. 

Two other processes produce peaks in the photon en- 
ergy spectrum close to the signal region: ISR produc- 
tion of the T(1S) and double radiative decays of the 
T(2S). The second transition in the processes T(2S) — > 
1X bj(lP),Xbj(lP) -► lT(lS), J = 0,1,2, produces 
peaks centered at 391, 423, and 442 MeV, respectively. 
These three peaks are merged due to photon energy reso- 
lution and the small Dopplcr broadening that arises from 
the motion of the Xw(l-P) m the cm. frame. We use 
the Xfcj(l-P) — * jY(lS) signal to validate estimates of 
signal efficiencies and determine the absolute photon en- 
ergy scale. Radiative production of the T(IS) via ini- 
tial state radiation (ISR), e + e _ — > JisrT(1S), leads to 



a peak near 550 MeV. The signal peak is better sepa- 
rated from the peaking background, with respect to the 
T(35) — > 7776 analysis [l|, primarily due to better abso- 
lute energy resolution at lower energy. 

Decays of the r/b via two gluons, expected to be its 
dominant decay mode, have high charged-particle mul- 
tiplicity. We select hadronic events by requiring four or 
more charged tracks in the event and that the ratio of the 
second to zeroth Fox- Wolfram moments [3] be less than 
0.98. 

Photon candidates are required to be isolated from all 
charged tracks. To ensure that their EMC shapes are 
consistent with an electromagnetic shower, the lateral 
moment [8[ is required to be less than 0.55. To ensure 
high reconstruction efficiency and good energy resolution, 
the signal photon candidate is required to lie in the cen- 
tral angular region of the EMC, —0.762 < cos(0 7i lab) < 
0.890, where # 7i lab is the angle between the photon and 
the beam axis in the laboratory frame. 

The correlation of the direction of the photon with the 
thrust axis 0] of the rjb is small, as there is no preferred 
direction in the decay of the spin-zero rjb and the mo- 
mentum of the rjb is small in the cm. frame. In contrast, 
there is a strong correlation between the photon direction 
and thrust axis in continuum events. The thrust axis is 
computed with all charged tracks and neutral calorime- 
ter clusters in the event, with the exception of the signal 
photon candidate. A requirement of Icos^tI < 0.8 is 
made to reduce continuum background, where 9t is the 
angle between the thrust axis and the momentum of the 
signal photon candidate. 

A principle source of background is photons from tt° 
decays. A signal photon candidate is rejected if in com- 
bination with another photon in the event it forms a tt° 
candidate whose mass is within 15 MeV/c 2 of the nominal 
ir° mass. To maintain high signal efficiency, we require 
the second photon of the ir° candidate to have an energy 
in the laboratory frame greater than 40 MeV. 

The final efficiency evaluated from simulated events is 
35.8%. 

The selection criteria were chosen by maximizing the 
ratio Ns/^/Nt, where Ns is the signal yield and Nt is 
the total yield of events in the signal region. The result of 
the optimization is insensitive to the exact definition of 
the signal region. A detailed Monte Carlo (MC) simula- 
tion [1 01 ] provides the signal sample for this optimization, 
while a small fraction (7%) of the T(2S) data is used to 
model the background. To avoid a potential bias, these 
data are excluded from the final fit of the photon en- 
ergy spectrum. The remaining T(2S) data used for the 
analysis have an integrated luminosity of 13 fb _1 , corre- 
sponding to (91.6 ± 0.9) million T(2S) events. 

To extract the 77b signal, a \ 2 fit of the E 1 spectrum 
is performed in the region 0.27 < E 1 < 0.80 GeV. The 
fit includes four components: nonpeaking background, 
Xb.j{^-P) -^7T(1S), jisrT(1S), and the rj b signal. 
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The nonpeaking background is 
an empirical probability density 



parametrized by 
function (PDF) 



xio 3 



where x = and A,C{ are deter- 



ii] 



^ ex p(Eti 

mined in the fit. 

Doppler-broadened Crystal Ball (CB) functions 
are used as phenomenological PDFs for the three 
Xb,/(1P) — > jT(1S) shapes. The CB function is a Gaus- 
sian modified to have an extended power-law tail on the 
low (left) side. The power law parameter describing the 
low-side tail of the CB function is common to all three 
of the Xbj(lP) peaks. The Doppler broadening of the 
X6j(lP) peaks is modeled by analytically convolving the 
CB functions with rectangular functions of half-width 
6.5, 5.5, and 4.9 MeV for the J = 0,1,2 states, respec- 
tively. These values are evaluated using the T(2S) and 
Xbj(lP) masses [H|. The resolution parameter of the 
Xbo(lP) PDF is fixed to that of the XbiO-P)- Due to 
its small yield and its position on the low side of the 
XbjO-P) peak, the exact width of the x&o(lP) is n °t 
crucial. The relative rates of the XbjO-P) components 
are fixed to values determined from a control sample of 

r(25) - ixbj(iP),xbj(iP) - jr(is),r(is) - 

events, and the relative peak positions from the world- 
averaged (PDG) values [121] . with a photon energy scale 
offset determined in the fit. 

The PDF of the peaking background from ISR T(1S) 
production is parametrized as a CB function with pa- 
rameters determined from simulated events. The ISR 
peak position is fixed to the value determined by the 
T(1S) and T(2S) masses [HI, minus the energy scale off- 
set shared with the x&j(1-P) peaks. 

The rjb signal PDF is a nonrelativistic Breit-Wigner 
function convolved with a CB function to account for the 
experimental _E 7 resolution. The CB parameters are de- 
termined from signal MC. Theoretical predictions for the 
r)b width, based on the expected ratio of the two-photon 
and two-gluon widths, range from 4 to 20 MeV [13] . Since 
the width of the r\b is not known, we have chosen a nom- 
inal value of 10 MeV, as in the T(3S) analysis. 

The free parameters in the fit are the normalizations 
of all fit components, all of the nonpeaking background 
PDF parameters, the 77b peak position, the energy scale 
offset, the Xbi(l-P) and x&2(lP) CB resolutions, and the 
transition point between the Gaussian and power law 
components of the Xbj(lP) CB functions. 

Figure 1 shows the photon energy spectrum and the 
fit result before (a) and after (b) subtraction of the non- 
peaking background. The x 2 per degree of freedom from 
the fit is 115/93. The line shapes of the three peaking 
components, %6j(lP), ISR T(15 l ), and the rjb signal, are 
clearly visible in the subtracted spectrum. The r]b signal 
yield is 12800 ± 3500 events, and the rjb peak energy is 
607.9^45 MeV. The observed signal width is consistent 
with being dominated by the resolution of 18 MeV. 

The ISR T(1S) yield can be estimated from data col- 
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FIG. 1: (a) Inclusive photon spectrum in the region 0.27 < 
E- t < 0.80 GeV. The fit is overlaid on the data points, (b) 
(color online) Inclusive photon spectrum after subtracting the 
nonpeaking background, with the PDFs for Xbj(lP) peak 
(light solid line), ISR T(1S) (dotted line), rjt signal (dashed 
line) and the sum of all three (dark solid line). 



lected below the T(AS) resonance. After correcting for 
the luminosity ratio, and the difference in ISR cross sec- 
tion and detection efficiencies at the two energies, we 
expect 16700 ± 700 ± 1200 ISR T(1S) events in the on- 
resonance T(2S) data sample. The consistency of the ob- 
served yield of the T(IS) component, 16800^^^ events, 
with the expected value provides an important valida- 
tion of the fitted background rate near the signal region. 
The yield and peak position of the 775 signal change by 
less than O.Ict when the ISR T(1S) yield is fixed to the 
expected value. 

We estimate the systematic uncertainty by varying 
the Breit-Wigner width in the rjb PDF to 5, 15, and 
20 MeV, varying the PDF parameters fixed in the fit 
by ±1 <t, using alternative smooth background shapes, 
varying the histogram binning between 1 and 15 MeV, 
incorporating a high-side tail to the Xbj(lP) peaks, and 
subtracting possible peaking background components. 
Smooth background PDF variations consist of using al- 
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ternative smooth background shapes that either incor- 
porate a 3rd-order polynomial in the exponential (i.e. 
c 4 = 0) or use a PDF of the form k{E 1 /E Q )~ Vl [1 + 
(£ 7 /i;o) 1/a ]~ (r2 ~ ri)Q - Other background shape varia- 
tions consisting of adding a term c§x 5 to the exponential 
of the smooth background function or adding a constant 
background PDF were found to change the fit negligi- 
bly. An additional high-side tail in the XbjO-P) peak 
may be produced by the coincidental overlap of photons 
from XbjO-P) decays with particles from the rest of the 
event or beam debris. We model this tail as a 90MeV 
wide Gaussian centered about each of the XbjO-P) peaks. 
Due to the large width of this component, it is indistin- 
guishable from the nonpeaking background, and its in- 
clusion does not improve the fit. We take the difference 
between the nominal fit and the fit including this tail 
as a systematic error. To evaluate the systematic due 
to the Xbo(lP) resolution, we perform a fit in which the 
XboO-P) resolution is fixed to that of the J = 1 state. 
To investigate the possible effect of peaking background 



from T(2S) 



,J])T(1S) events, these contributions 



are subtracted prior to fitting, assuming the measured 
value and 90% CL upper limits for the branching frac- 
tion of the 77 and tt° [12 1 transitions, respectively, giving a 
variation of —71 (+651) events for the 77 (tt°) transition. 
Photons from the transition T(2S) — > 7r°7r T(15), which 
produce a smoothly varying background below 400 MeV, 
are absorbed into the smooth background PDF and do 
not require a separate systematic error. 

Including systematic uncertainties, the signal yield is 



12800 ± 3500 



+3500 
-3100 



events. The largest contributions to 



the systematic error on the yield are from the ?7b width 
variation (l^oo events) and the background shape vari- 
ation (^2700 events). 

The photon energy scale is corrected with the fitted en- 
ergy scale offset of 1.4±0.2±0.7 MeV determined from 
the Xbj(lP) and ISR peaks. The systematic error is half 
of the shift added in quadrature with the PDG errors on 
the Xbj(lP) masses (0.4 MeV). The ISR peak contributes 
negligibly to the determination of the offset, due to its 
small yield. After including an additional systematic un- 
certainty of 1.7 MeV from the fit variations described 
above, we obtain a value of E 1 = 609.3±f% ± 1.9 MeV 
for the rjb signal. 

To confirm that this state is identical to the state 
observed in the Y(3S) — > 7775 analysis [l| we calcu- 
late the significance of the signal with the signal peak 
fixed to 614.3 MeV, the value expected for an rjb mass of 
9388.9 MeV. The r\b signal significance is estimated as 
iV^ 2 (no signal) — x 2 (fixed mass), where % 2 (fixed mass) 
is the x 2 °f the fit with the % signal included and 
X 2 (no signal) is the x 2 of the fit with the r\b PDF re- 
moved. The statistical significance estimated in this way 
is 3.7 standard deviations. The significance of the sig- 
nal, including systematics, is estimated by making the 



variations discussed above. Additional cross-checks are 
performed by changing the lower (upper) limit of the fit 
range to 250 MeV (850 MeV) and varying the selection 
on I cos(#t)|. In all fits, the significance lies between 3.0 
and 4.3 standard deviations. 

The rjb mass derived from the E 1 signal is M(rjb) = 
9394.2t^ ± 2.0 MeV/c 2 . Using the PDG value of 
9460.3 ± 0.3 MeV/c 2 for the T(15) mass, we determine 
the r(lS)-r) b mass splitting to be 66.ll 4 ;^±2.0 MeV/c 2 . 

For the measurement of the branching fraction, we 
have an additional source of uncertainty resulting from 
the signal selection efficiency. The systematic uncer- 
tainty on the photon detection efficiency is 1.8%. We 
estimate the uncertainty on the hadronic selection ef- 
ficiency (4.9%) by comparing data and MC efficiencies 
of the selection on hadronic T(1S) events. The uncer- 
tainty in photon quality selection efficiency (0.5%) is es- 
timated from tt° decays in data and MC. The difference 
between the efficiency in MC and the efficiency for a 
flat distribution (0.6%) is used as the uncertainty on the 
I cos 9t I selection. We determine the uncertainty for the 
7T° selection (4.1%) by comparing the efficiency-corrected 
XbjO-P) yield with and without the 7r° veto. The total 
systematic error on the selection efficiency is 6.7%. The 
uncertainty on the number of T(2S) events is 0.9%. In- 
corporating these systematic uncertainties, we determine 
the branching fraction of the decay T(2S) — ► jrjb to be 
(3.9±l.lt^)xl0- 4 . 

In the T(3S) analysis [lj, we estimated the systematic 
uncertainty on the signal efficiency using Xbj(2P) de- 
cays, incurring a large error (22%) due to the uncertain- 
ties in the Xbj(2P) branching fractions. The uncertainty 
in T(3S) — * 7776 efficiency obtained using the procedure 
described above is 5.5%, resulting in a final branching 
fraction of B[Y(3S) -> jrj b ] = (4.8 ± 0.5 ± 0.6) x 10~ 4 . 
This value supersedes our previous result, which differs 
only in having a systematic uncertainty two times larger. 

Using the results given above, we determine a branch- 
ing fraction ratio of B[T{2S) -> -yr) b ]/B[r(3S) -> Tib] = 
0.82 ± 0.24t°; 2 °. The systematic uncertainties due to 
selection efficiency and the unknown r\b width partially 
cancel in the ratio. Our measurement is consistent with 
some of the theoretical estimates of this ratio of magnetic 
dipole transitions to the 77/,, 0.3 — 0.7 2f, while the abso- 
lute transition rates are not well-predicted by theoretical 
models. 

In conclusion, we have obtained evidence, with a signif- 
icance of 3.0 standard deviations, for the radiative decay 
of the T(2S) to a narrow state with a mass slightly less 
than that of the T(1S). The ratio of the radiative pro- 
duction rates for this state at the T(2S) and T(3S) reso- 
nances is consistent with that expected of the 1%. Under 
this interpretation, the mass of the 77^ is 9394. 2l4g ± 2.0 
MeV/c 2 , which corresponds to a mass splitting between 
the r(lS) and the 776 of 66.1±|;| ± 2.0 MeV/c 2 , consis- 
tent with the value from the T(3S) analysis. The aver- 
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age of the two results is M(r) b ) = 9390.8 ± 3.2 MeV/c 2 . 
This value of the r)b mass is consistent with a recent un- 
quenched lattice prediction [j| but more than two stan- 
dard deviations away from the mass predicted by ap- 
proaches based on perturbative QCD 
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